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Abstract 

The higher-order superintegrability of systems separable in polar coordinates is 
studied using an approch that was previously applied for the study of the superin- 
tegrability of a generalized Smorodinsky-Winternitz system. The idea is that the 
additional constant of motion can be factorized as the product of powers of two 
particular rather simple complex functions (here denoted by M and N). This tech- 
nique leads to a proof of the superintegrability of the TremblayTurbinerWinternitz 
system and to the explicit expression of the constants of motion. A second family 
(related with the first one) of superintegrable systems is also studied. 
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1 Introduction 



A superintegrable system is a system that is integrable (in the sense of Liouville- Arnold) 
and that, in addition to this, possesses more constants of motion than degrees of freedom. 
Most of known superintegrable system are systems which admit separability in two (or 
more) different coordinate systems; for example, the so-called Smorodinsky-Winternitz 
(SW) potential 

t r 12/2, 2\ , ^2 &3 

V sw = -u (x +y ) + -2 +-2 , 

jj y 

is separable in Cartesian and polar coordinates [I] [3]. A very important point is that the 
three constants of motion are quadratic in the velocities (momenta). In fact, if we call su- 
perseparable a system that admits Hamilton- Jacobi separation of variables (Schrodinger 
in the quantum case) in more than one coordinate system, then quadratic superintegrabil- 
ity (i.e., superintegrability with linear or quadratic constants of motion) can be considered 
property arising from superseparability. 

On the other hand, it is known the existence of some systems endowed with 'higher-order 
superintegrability', that is, with integrals of motion which are polynomials in the momenta 
of order higher than two (see Refs [I]- [21] for some papers published in these last years). 
Some of them are not separable in any system of coordinates (as the Calogero-Moser 
system j22]-[25j), but others are separable in only one particular system of coordinates. 
In this second case, the system possesses two quadratic integrals (arising from separability) 
and a third integral of higher order. 

We will focus our attention on two higher order superintegrable systems related with 
the SW system. 



(i) A natural generalization of V sw is given by the following potential 

t t / \ 1 9/ 9 9 9 9\ 1 ^2 

V(n x , n y ) = -u {n x x + n y y ) + — + — , 

that contains a more general harmonic oscillator with anisotropy [25]— [29]. When 
n x — n y — 1 it becomes the SW system, but in the general case it is only separable 
in Cartesian coordinates; so, the third integral is a polynomial in the momenta of 
higher order than 2. 

(ii) The following potential 

V ttw (r, 0) = | w V + ( — + . z ) , 
1 2r z \cos z (k(p) sm (/c0)/ 

has been firstly studied by Tremblay, Turbiner, and Winternitz [30] [31] , and then 
by other authors [32]-[36j. When k — 1 it becomes the SW system, but in the 
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general k ^ 1 case it is only separable in polar coordinates; therefore, the third 
integral must be a polynomial in the momenta of higher order than two. 

The superintegrability of the potential V(n x ,n y ) has been studied using different ap- 
proaches [26], [27J, [28j EH]. We point out that it was proved in [29] the existence of a 
complex factorization for the higher-order constant of the motion. This approach was 
in fact a generalization of a method previously used (37J EE] for the superintegrability of 
the harmonic oscillator with commensurable frequencies (the ratio is a rational number). 
Concerning the TTW system it has been studied at both the classical and the quantum 
levels. We quote here some different approaches: (i) Proof that the quantum system is 
exactly solvable and that it is endowed with certain properties that represent a very strong 
indication of superinegrability [30J; (ii) Proof that all bounded classical trajectories are 
closed [31]; (hi) Techniques related with supersymmetric quantum mechanics [32]-[33]; 
(iv) Hamilton- Jacobi formalism [31]; (hi) As a particular case arising in the general study 
of the Schrodinger equation for quantum systems that admit a separation of variables in 
some orthogonal coordinate system. [35] • We also mention that it has been the starting 
point for the study [22] , [IS] , [2S] , of quantum algebraic properties related with intertwin- 
ing operators. Very recently, Gonera has proved (in a very modified version of [34]) the 
superintegrability of this system making use of the formalism of the angle-action variables 

PI- 

We also mention that Post and Winternitz have studied in [40] the quantum superinte- 
grability of a family of Hamiltonians with a deformed Kepler- Coulomb potential, and also 
depending of a parameter k, by relating this new family with the TTW system using the 
so-called coupling constant metamorphosis (or Stackel transform) [41] as an approach. 

The purpose of this paper is twofold: to study the existence of higher-order superinte- 
grable systems that are separable in polar coordinates using as an approach the method 
applied for the study of the superintegrability of V(n x ,n y ) and to study the superinte- 
grability of the Tremblay-Turbiner- Winternitz system. 

The structure of this article is as follows: in Sec. 2, we present the main characteristics 
of the method in the Cartesian case; we first consider the linear harmonic oscillator and 
then we recall the proof of the superintegrability of the nonlinear system V(n x ,n y ). In 
Sec. 3, we study the higher-order superintegrability of a family of potentials separable 
in polar coordinates and we prove the superintegrability of the TTW system. In Sec. 4, 
we obtain a new family of superintegrable potentials obtained from the previous one by 
making use of a transformation (we also study the relation of this family with the TTW 
system). Finally, in Sec. 5, we make some comments and present some open questions. 
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2 Complex factorization and superintegrability : Sepa- 
rable systems in Cartesian coordinates 

In this section we recall the method of the complex factorization for two oscillators sepa- 
rable in Cartesian coordinates. 

2.1 Superintegrability of the linear oscillator 

The two-dimensional harmonic oscillator 

L HO = \{v 2 x + vl)- l -{u l 2 x 2 + uj 2 W) (1) 

has the two one-degree of freedom energies, ii = E x and I 2 = E y , as fundamental integrals. 
The superintegrability of the rational case, u)\ = n x co , U2 = n y uo, with integers n x ,n y , 
can be proved by making use of a complex formalism [37J [38]. Let A4, i = x,y, be the 
following two complex functions 

A x = v x + i n x u x , A y = v y + i n y u y , 

then we have the following time-evolution 

d A A d A A 

—A x = m x u} A x , —A y = in y oj A y . 

Thus, the functions Ay defined as 

A l3 = {A^{A]) n % i,j = x,y, 

are constants of motion. The two real functions lA^a-l 2 and (A^l 2 are proportional to the 
energies E x and E y and concerning A xy , since it is a complex function, it determines not 
one but two real first integrals, lm(A xy ) and Ke(A xy ). So, we have obtained four integrals 
but, since the system is two-dimensional, only three of them can be independent. We can 
choose ^3 = lm(A xy ) as the third fundamental constants of motion (the other constant 
I4 = Ke(A xy ) can be expressed as a function of E x , E y , and Im^A^)). 



2.2 Superintegrability of a generalizated version of the SW sys- 
tem 

A natural generalization of the SW potential [TJ [2] is given by 

V{n x , n y ) = \u Q \nlx* + n\y>) + A + h- , ( 2 ) 
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representing an harmonic oscillator with rational ratio of frequencies and inversely quadratic 
nonlinearities. It is separable but not multiple separable and because of this the addi- 
tional constants of the motion are of higher order in the momenta. It has been proved in 
[21] that it also admits a complex factorization that was obtained as a deformation of the 
quadratic version of the factorization of the integrals of motion of the linear oscillator. 

The proof is as follows. Let us now denote by Bi, i — x,y, the following two complex 
functions 



B x = (v 2 x - n 2 x u 2 x 2 + + 2 i n x u xv x , B y = (v 2 - n 2 y u 2 y 2 + -^j + 2 i n y u yv y . 



then we have the following time-evolution 

— B x = 2i n x u) B x , —B y = 2i n y uo B y . 

Thus, the complex functions B^ defined as 

B tj = (Bip (B*p , i,j = x,y, 

are constants of the motion 

j t B xy = B^B^ [n y B x B; + n x B x £*) = . 

The moduli of B x and B y are related with the energies 

|^| 2 = A(E 2 X - hnlu 2 ) , \B y \ 2 = A(E 2 - k 2 n 2 y ul) . 
and, as in the linear case, B xy can be considered as coupling the two degrees of freedom. 



3 Complex factorization and superintegrability : Sepa- 
rable systems in polar coordinates 

The two systems studied in Sees. (12. ip and (I2.2p were separable in Cartesian coordinates 
and the two functions, A{ in the linear case and Bi in the nonlinear case, were dependent 
of only one coordinate. Now, in this section, we study the existence of a complex factor- 
ization but for potentials separable in polar coordinates (r, 0). The main difference is that 
now we present the inverse approach; we start with a general separable potential V(r, <fi) 
and then we determine the conditions to be satisfied by the potential to be superintegrable 
with the additional constant of motion given by the product of two complex functions. 

We will make use of the Hamiltonian formalism, so, in what follows, the time derivative 
d/dt means Poisson bracket with the function H. 
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Let us consider the following Hamiltonian 

tt 1 / 2 , Pl\ ,1 2 2,1 F{4>) 

H = 2^ + ^+2^ T + 2— ■ ( 3 ) 

It is separable in polar coordinates and it is therefore endowed with the following two 
constants of the motion 

2 , Pi , , . 2 2 , F(<t>) 



h = p; + ^ + woV + 

J 2 = pl + F(cf>) 



Any property to be satisfied by a Hamiltonian as (J3j) must be true for the case of the 
single oscillator (without additional ^-dependent term). So we first consider the case of 
a central potential without the function F. 

Step 1. Let us first consider a central potential U(r). Then if we denote by M r and Mj 
the functions 

2 



M r = - PrP4> , Mi=p 1 + tf(r) 



then we have 



(i) The time-derivative of Mj is given by 



r ■ 



(ii) The time-derivative of M r is given by 

|* = -2(i W )*i, 
if and only if the potential C/ is solution of the following equation 

rU' -2U = 0. 

The harmonic oscillator satisfies the properties (i) and (ii) with the factor Ao given 
by Ao = (p^/r 2 ). It seems natural to assume that, when considering the more general 
Hamiltonian ([3]), then the momentum p^ in M T and Ao will probably change to \fJ2~. 

Step 2. Consider now the non-linear potential 



T/Y A\ 1 2 2, lF{4>) 
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Then if we denote by M r and M, the functions 



2 2 

M r = - Pr V J 2 , Mi = J 1 - — J 2 . 

the time-evolution of them is given by 

^-M r = 2(-X)M l , j t M t = 2XM r , \ = ^^J 2 . 

Step 3. Consider a potential separable in polar coordinates. 



V (r,(f)) = U(r) + l 

and denote by N r the following function 

N r = z + J2 cos(A;0) , 
where z denotes a real parameter. Then we have 

(i) The time-derivative of N r is given by 



L N r = (-A) k Ni , N i = \[J 2 v<t, sin(/c0) , A = ^ 



(ii) The time-derivative of Ni is given by 

^-Ni = \kN r , 
at 

if and only if the function F is solution of the following equation 

sin(A;0) F' + 2k cos(A;0) F + 2kz = . 



This is a linear first order equation for F and the general solution is given by 

F = ka + 2z ( cos ( k ^ ) 
sin 2 (/c0) Vsin 2 (A;0)/ ' 

Notice that F is a linear combination of two functions because the presence of the param- 
eter z. If z = then F reduces to the solution of the homogeneous version of the linear 
equation. 

The steps 1 and 2 provide a property that singles out the harmonic oscillator in the set of 
all central potentials; on the other hand, the step 3 states the existence of an appropriate 
function iV for any function U(r) and certain very particular values of the function F(4>). 
The superposition of these two properties leads to the following proposition. 
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Proposition 1 Consider the non-linear potential 

\r( a\ 1 22, i-W) j? K fcos{k(f))\ 

where k a and kb are arbitrary constants. Let M and N be the complex functions M = 
M r + i Mi and N = N r + i TV, with real and imaginary parts M a and N a , a — r,i, be 
defined as 

M r = ^p r ^/T 2 , M i = J 1 -^J 2 , 
k 

N r = -^ + J 2 cos(/c0) , Ni = \[7 2 p (i> sin(A;0) . 
Then, the complex junction K defined as 

K = M k (N*f 
is a (complex) constant of the motion. 

Proof: First, let us comment that if k is an integer (A; = n) then V\ is well defined as 
a function on the plain and the domain (configuration space of the system) is the union 
of 2n wedges that cover the plain (the singularities of V\ are in the borderlines). If k is 
a rational number then the configuration space must be restricted to only one wedge of 
angle smaller than 2tt. 

The time-evolution of the functions M and N is given by 

4-M = i2AM, ^-N = ik\N, A = 4 
dt dt r< 



Thus we have 



^-K = M^-^N* ( kMN* + 2MN*) = 0. 
dt V / 



d_ 

dt 

□ 

Next, we give the expressions of the function F\ for the four first integer values of k: 

ka kfc x 



Fi(k 


= 1) 


F^k 


= 2) 


F,(k 


= 3) 


F^k 


= 4) 



y 2 y 2 ^x 2 + y 2 
k a kb k a + kb 

Ax 2 Ay 2 

0,2 _! y2) 2,2 { k ^ 2 + y 2 ) 2 + ^ - 3y 2 )*V^Tf) 



8 



At this point we recall that the potential of Tremblay-Turbiner-Winternitz [30] [36] is 
given by 

V ttw (r, 4>) = \ Wo 2 r 2 + 7r~2~ f — + • Jms ) • 
1 2r A \cos z (kcp) sin (hep)/ 

It can be proved that the function F\ reduces to the ^-dependent part of the TTW 
potential when (i) the coefficients k a and kb are writen as k a = 2(a + /3) and k = 2((3 — a), 
and (ii) the angular coefficient k is changed to 2k; that is, 

2(a + ff) (3 , / cos(2/c0U a 

sm\2k4>) [P ' Vsin 2 (2A;0)/ cos 2 (£;0) sin 2 (fc0) ' 

and consequently 

Vi(r, cp, 2k) = V ttw {r,(f), k) . 

In fact, we must point out that the expression obtained for K in proposition [1] is quite 
similar (introducing the appropriate changes of coefficients and notatios) to the expression 
obtained by Gonera in [39J making use of the action-angle formalism. 

Summarizing, the above proposition [1] states that the potential V\, that is a separable 
potential, it is also superintegrable with a third constant of the motion that arises as the 
real or the imaginary part of a complex function K obtained from the functions M and 
N. Moreover, this potential is equivalent to the TTV potential; therefore, proposition [1] 
also proves the superintegrability of the TTW system. 



4 A second family of superintegrable systems 

Suppose that we interchange the cos(k<f)) and the sin(k<f)) functions in the expression of 
the function N. In this case we then we arrive (in the step 3) to the equation 

cos{k<p) F' + 2k sin(A;0) F + 2kz = , 

wich has as general solution the function F2 given by 



cos 2 (kcp) V cos 2 (kcp) J 



2z) . 



Consequently, we can also state that the potential 

t/ i sk\ 1 22, 1-^(0) „ k a /sin(kcp)\ 

is superintegrable with a third (complex) constant of the motion given by 



K = M k (Np 2 



2) 5 
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where M is defined as in proposition 1 and N 2 is now defined asiV 2 = N 2r + % N 2i with 
N 2r = — + ^ 2 cos (^0) , N 2i = \fhv<i> sin(/c0) . 
Next, we give the expressions of the function F 2 for the four first integer values of k: 
F 2 (k = 1) 

F 2 (k = 2) 
F 2 (k = 3) 
F 2 (k = 4) 

Finally, we consider the relation of F 2 with F\ and with TTW. 

(i) The interchange the cos(/c0) and the sin(/c0) functions can also be seen as a rotation 
of angle n/(2k). That is, 

F 2 (0^0 + vr/(2A;)) = F 1 . 

The kinetic term in H is rotationally invariant; so any rotation (arbitrary value 
of the angle) transforms the potential V\ in a new potential that is therefore also 
integrable (the new constants are the rotated versions of the original constants). 
From this viewpoint, the potential V 2 can be seen not just as a new potential but as 
a rotated version of V\. Nevertheless the transformation V\ — > V 2 is rather peculiar 
and interesting since it is a different rotation for every value of k and is such that 
it preserves the high degree of simplicity of the potential. We also note that the 
domain of V 2 is different to the domain of V\. 

(ii) We have proved that 

V M = V ttM ,2) = I W0 V + ±- 2 (— ^ + T -J^) . 

Then if we replace cos(/c0) by sin(/c0), we obtain the following TTW-related poten- 
tial 

Vttwir, <P,k) = -ujq r + 7r ^[ : 7TT\ + i ■ 7T~T\ ■ 

2 2H VI + sm(£;0) 1 — sm(fc0)/ 

So that we arrive at the following relation 

\/ 2 (r,0, fc) = V ttw (r,(f), k) . 
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x 2 + y 2 2xy 



\/x 2 + y 2 



(x 2 — I/ 2 ) 2 



X 



(a; 2 — 2/2)2 

2 _ ^, 2VU2 (k a (x 2 + y 2 ) 2 + k b (3x 2 - y 2 )y^T7) 



3y 2 ) 2 x 2 
1 



(x 4 - Qx 2 y 2 + y 4 ) 2 



k a (x 2 + y 2 f + Ak b {x A - y 4 )xy 



5 Comments and open questions 



We have studied the existence of systems that are separable in polar coordinates and 
admit an additional constant of motion of higher order in the momenta. The main idea is 
that this additional constant can be factorized as the product of powers of two particular 
rather simple complex functions (here denoted by M and N). An important point is that 
all the potentials obtained by this method are related to the harmonic oscillator. As is 
well known, the other fundamental superintegrable central potential is the potential of the 
Kepler problem (this is just the theorem of Bertrand). The Kepler problem admits two 
noncentral superintegrable potentials separable in two different systems (this was proved 
in p]). Moreover, the existence of deformations with higher order superintegrability has 
been proven recently [ID] . It would be convenient try to apply the procedures presented in 
this study to the Kepler problem. This question, as well as some other related problems 
(as the properties of the quantum version of the functions M and N), is an open question 
to be studied. 
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